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ABSTRACT 

We present an X-ray image of the BL Lacertae object OJ287 revealing a long jet, curved 
by 55° and extending 20", or 90 kpc from the nucleus. This de-projects to > 1 Mpc based 
on the viewing angle on parsec scales. Radio emission follows the general X-ray morphology 
but extends even farther from the nucleus. The upper limit to the isotropic radio luminosity, 
~ 2 x 10 24 W Hz -1 , places the source in the Fanaroff- Riley 1 (FR 1) class, as expected for 
BL Lac objects. The spectral energy distribution indicates that the extended X-ray emission 
is from inverse Compton scattering of cosmic microwave background photons. In this case, the 
derived magnetic field is B « 5 /iG, the minimum electron energy is 7-40m c c 2 , and the Doppler 
factor is S ~ 8 in a knot 8" from the nucleus. The minimum total kinetic power of the jet is 
1-2 xlO 45 erg s _1 . Upstream of the bend, the width of the X-ray emission in the jet is about half 
the projected distance from the nucleus. This implies that the highly relativistic bulk motion is 
not limited to an extremely thin spine, as has been proposed previously for FR 1 sources. The 
bending of the jet, the deceleration of the flow from parsec to kiloparsec scales, and the knotty 
structure can all be caused by standing shocks inclined by ~ 7° to the jet axis. Moving shocks 
resulting from major changes in the flow properties can also reproduce the knotty structure, but 
such a model does not explain as many of the observational details. 
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1. Introduction 

According to the unified scheme of extragalac- 
tic radio sources, BL Lacertae (BL Lac) objects 
repre sent Fanaroff- Riley 1 (FR l uFanaroff fe Rilev 
19741 ) radio galaxies in which o ne of the jets points 
within ~ 15° of the line of sigh t (|Antonucci fe Ulvestad 
19851: iPadovani fc Urrvlll990l) . Based on a care- 
ful analysis of the radio jets of 38 radio sources, 
Laing et al. (Il999l) inferred that the bulk flow of 
a typical FR 1 jet decelerates from a relativistic 
speed to less than ~ 0.1c within a distance of two 
to tens of kpc from the nucleus. The inferred de- 
celeration agrees with the expectations of a model 
in which t he jet entrains material fro m its sur- 
roundings (|Bicknel]|[l994l lLaing|[i996l) . The dis- 
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tance from the nucleus over which the deceleration 
occurs is then expected to increase with the kinetic 
power of the jet. It should therefore be greater 
at higher values of both the initial bulk Lorentz 
factor and the luminosity of extended emission. 
Laing et al. (1999) also found that the data are 
consistent with a model in which the velocity of 
the flow decreases from the axis to the edge of the 
jet, either gradually or in a more distinct manner , 

(2000) 



Chiaberge et al 



as a spine plus sheath, 
came to the same conclusion based on a compar- 
ison of BL Lac and FR I core spectral energy 
distributions. Some jet launching models include 
an ultra-fast spine produced by the ergosphere of 
a rotating black hole, with the slowe r sheath ema- 
nating from the accretion disk ( e.g., Punslv [l996: 



McKmnevlkoOfjl) . Alternatively, iMeierl (|2003[ ) pro- 



posed that the spine is generated by the central 
engine, while the sheath arises from velocity shear 
downstream. 

The proposed deceleration of jets could rec- 
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oncile gross discrepancies between inferred bulk 
Loren tz factors T of the BL Lac and FR I popula- 
tions. Urrv. Padovani. fe Stickell ( 1991 ) concluded 
that the properties of BL Lac objects require the 
values of T of the FR I parent population to range 
up to ~ 30. There is independent support for 
this conclusion: measureme nts of proper motions 
of 7-ray bright BL Lacs bv Ijorstad et al. ( 2001 ) 
indicate apparent speeds (and therefore minimum 
values of T) as high as ~ 40c. 

The extended X-ray emission from low-power 
FR 1 radio galaxy jets is consistent with syn- 
chrotron radiation within the inner 10 kpc by 
electrons with energies in the 10-100 TeV range 
(IWorrall. Birkinshaw. fc Hardcastlel200ll ; lKataoka 



as is the case for some quasars (jDermer fc Atovan 



2004). Furthermore, if the jet possesses a spine- 



2005). Those that are more luminous, highly rel- 
ativistic, and pointing close to the line of sight 
- characteristics of radio-bright BL Lac ob- 
jects — can also produce observable X-rays from 
inverse Compton scatter ing off the mic r owave 
background (IC/CMB; iTavecchio et all [2 000; 



Celotti. Ghisellini. fc Chiaberge); 2001; Sambruna et al 



200 ll ) on scales of tens or hundreds of kiloparsecs. 
The total energy requirements are higher in the 
IC / CMB case owing to the radiative inefficienc y 
of low-energy electrons ([Dermer fc Atovan|[200l . 
This has led to a controversy over the process by 
which X-ray emission is produced in the jets of 
quasars and BL Lac objects with highly relativis- 
tic jets. 

Until now, observations that provide informa- 
tion on deceleration, spine-sheath structure, and 
X-ray emission have been limited to FR 1 ra- 
dio galaxies that are not BL Lac objec ts (e.g., 
Worrall. Birkinshaw. fc Hardcastlel 2001 ) , three 
BL Lac objec ts with mixed FR 1 and FR 2 
morphologies (jSambruna et al.l 120071 120081 ). and 
blazar-class quasars with powerful, high ly rel- 



ativistic jets (e.g.. i Jorstad fc Marscherl l2004t 
Harris fc Krawczvnski 20061 ). The extended X- 
ray emission of low-power, high-r BL Lac objects 
has been less well studied. Such jets should ra- 
diate at X-ray energies via the IC/CMB and/or 
synchrotron process, with the emission gradually 
becoming weaker with distance from the nucleus 
as the jet decelerates. Because of the low power of 
the jet as determined by the luminosity of the ex- 
tended emission, the IC / CMB process can be elim- 
inated as a major contributor to the X-ray emis- 
sion if the total energy requirement is too high, 



sheath structure, the IC/CMB emission should 
arise only in the narrow, highly relativistic spine. 

In order to explore the low-power, high-r case, 
we have carried out multi-waveband imaging of 
the jet of the BL Lac object OJ287 with the Chan- 
dra X-ray Observatory (CXO), Hubble Space Tele- 
scope (HST), Spitzer Space Telescope (SST), and 
Very Large Array (VLA). We report here the de- 
tection of X-ray emission along a length of the 
jet extending ~ 90 kpc from nucleus in projec- 
tion, and more than 1 Mpc after de-projection. 

The X-ray emission, which has a somewhat knotty 

fc Stawjefcra cture. is resolved in the transverse direction 
and fades gradually with distance from the nu- 
cleus, as does the radio emission. The radio and 
X-ray fluxes, along with upper limits at infrared 
and optical wavelengths, provide a spectral energy 
distribution that strongly constrains the emission 
model. 

We discuss the observations in §2 and analyze 
the results in §3. §4 contains our interpretation 
of the observations and discusses the implications 
for extended jets, while §5 presents our main con- 
clusions. In this study, we adopt the current stan- 
dard flat-spacetime cosmology, with Hubble con- 
stant #0=71 km s" 1 Mpc" 1 , tt M = 0.27, and 
n A = 0.73. At a redshift z = 0.306, OJ287 has 
a luminosity distance dg_ = 1.58 Gpc, and 1" cor- 
responds to a projected distance of 4.48 kpc in the 
rest frame of the host galaxy of OJ287. 

2. Observations and Data Analysis 

Table Q] lists the main information about the 
CXO, HST, SST, and VLA observations of OJ287. 
The VLA data include both a new observation at 
a wavelength of 2 cm and archival data at 6 and 
21 cm. 

2.1. X-Ray Observations 

We observed OJ287 using the back-illuminated 
S3 chip of the Advanced Camera for Imaging and 
Spectroscopy (ACIS) on CXO and a 1/8 subar- 
ray to reduce the frame time to 0.4 s. The data 
analysis, which utilized version 3.3.0 of the CIAO 
software and version 3.2.4 of CALDB, followed the 



procedures detailed in lJorstad fc Marscherl (|2004l ) . 
We generated a new level 2 event file with a pixel 
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size of O'.'l and convolved the image with a 0'.'5 
FWHM circular Gaussian kernel corresponding to 
the angular resolution of the ACIS. The resulting 
image (Fig. [I} shows a prominent jet that extends 
out to ~ 20" from the nucleus. The direction of 
the jet changes by 55° from the straight inner sec- 
tion within 8" of the core to the end of the X-ray 
emission. The jet has a knotty structure, which 
we separate into five regions, Jl — J5. All knots 
are detected at a level of 10cr or higher relative to 
the background. 

We fit the spectral data for the jet in the pho- 
ton energy range 0.2-6 keV by a power-law model 
with fi xed Galactic absorption N g — 3.02 x 10 20 



cm (jDickev &: Lockmanl 119901) . (Allowing the 



column density to vary yields a value that is not 
significantly different from this.) The calculated 
parameters of the X-ray jet are given in Tabled 

2.2. Optical Observations 

Over two orbits of HST, we obtained 12 expo- 
sures of 230 sec each of OJ287 with the WFPC2 
PC-chip using the broad filter F606W in a 4-point 
dithering mode. The images were combined with 
task crrej within the STSDAS software package. 
We generated the point-spread function (PSF) 
with the TinyTim software and subtracted it from 
the combined image. The deconvolved image is 
displayed in Figure [5] We did not detect signifi- 
cant flux over the sky level along the position of 
the X- ray and radio jet. lYannv. Jannuzi. fc Impev 
(|1997j ) also did not detect the jet at 814 nm with 
the WFPC2, and found that emission from nebu- 
losity associated with the host galaxy is insignifi- 
cant beyond 4" from the nucleus. 

2.3. Infrared Observations 

We observed OJ287 with IRAC on SST at 3.6 
and 5.8 /im with a single frame exposure of 
12 s in a dithering mode (36 Positions Reuleaux). 
We used the Level 2 Post-BCD Pipeline images, 
which we checked for artifacts. We created the 
PSF at each infrared (IR) wavelength using im- 
age stars, then subtracted it from the image with 
the task apex_qa_lframe.pl of version 18.3.1 of the 
MOPEX/APEX software package. Figure [3] dis- 
plays the deconvolved image at 3.6 /xm, which re- 
veals some enhanced brightness in the jet direc- 
tion, especially in the regions of knots J2 and 



J4. However, the fluxes are at l-2er detection 
level owing to possible confusion with unrelated 
IR sources, hence we cannot ascertain that the 
emission is physically associated with the jet. We 
therefore use the measured fluxes (plus the 2-a un- 
certainties) to estimate the IR upper limits for the 
jet regions. 

2.4. Radio Observations 

We processed the VLA data with the Astro- 
nomical Image Processing System (AIPS) software 
provided by the National Radio Astronomy Ob- 
servatory (NRAO), following the standard proce- 
dures outlined in the AIPS Cookbook available at 
www.nrao.edu. (The observations at 15 GHz in- 
cluded only 15 of the usual 27 antennas, since 
they occurred during the changeover to the new 
EVLA.) Each observation employed two 50 MHz 
bands, centered on 1.3851, 1.4649 GHz (21 cm 
wavelength), 4.5351, 4.8851 GHz (6 cm), and 
14.915, 14.965 GHz (2 cm), with mean frequen- 
cies listed in Table _Q We performed the i maging 
with the Difmap software ( Shepherd~lll997l ) . again 
following standard procedures. Figure 0] presents 
all three radio images, while Table [2] gives the pa- 
rameters of the radio emission. The image at 1.425 
GHz is very similar to the full-resolu tion version 
presented bv lPerlman k, Stockd ( 1994 ). which was 
based on the same data. 

3. Observational Results 

Figure [5] displays the entire jet at both radio 
and X-ray frequencies. The radio and X-ray emis- 
sion is generally coincident within 12" of the nu- 
cleus, beyond which the radio emission falls below 
the detection limit. An exception is knot J3, sit- 
uated along the southern rim of the jet with cen- 
troid 7'.'5 from the nucleus, where we detect no 
radio emission above the noise level. Because of 
the low radio flux, we cannot reliably determine 
the flux density separately at 4.7 or 14.9 GHz for 
each knot denoted in Figure [1] We therefore mea- 
sure these quantities only for the entire extended 
jet from knot Jl to J5, with the values given in 
Tabled 
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3.1. Comparison of X-ray and Radio Emis- 
sion 

Figure [5] marks the jet axis (which traces the 
mid-point of the detected X-ray emission trans- 
verse to the axis) and three crosscuts that pass 
through knots J2, J3, and J4. In Figure [7] we 
plot both the X-ray and radio intensity profiles 
along the jet axis, and only the transverse X-ray 
profile along the crosscuts. The path along the 
axis misses knot J3, which lies along the south- 
ern boundary of the jet, but crosscut b includes 
the knot. Knots Jl, J2, and J4 are clearly seen in 
the axial profile, while crosscut a passes through 
knot J2 and c includes J4. The maximum X-ray 
intensity lies farther from the nucleus than does 
the radio peak in both Jl and J2 by 0'.'9 and 0"5, 
respectively, while the peaks at the two wavebands 
coincide in knot J4. In each case, however, the up- 
stream boundaries of the radio and X-ray emission 
coincide. 

The displacement of the intensity peaks of 
knots J2-J4 from the axis (defined as the centroid 
of the detectable X-ray emission across the jet) 
seen in Figures [T] and [7] is striking. On the other 
hand, the intensity along the axis is too strong to 
consider the jet to be edge-brightened. 

The spectral index of the radio emission from 
the jet, integrated over the region containing de- 
tectable X-ray flux, is a(radio)= 0.8 ± 0.1. This 
differs somewhat from the X-ray value, a (X-ray) = 
0.61 ±0.06. The discrepancy could be due to a gra- 
dient in X-ray spectral index, which increases with 
distance from the nucleus from knots J2 to J4. 
We note that the overall radio spectral index falls 
within the 1-a range of the X-ray spectral index 
for all knots except J2. Given this, and because 
we can only determine the overall radio spectral 
index, we will ignore the small apparent radio/X- 
ray discrepancy in spectral slope and adopt a value 
a = 0.7 in our analysis. 

The radio emission on milliarcsecond scales con- 
sists of a bright "core" with multiple knots - 
some with apparent superluminal motion and oth- 
ers that are subluminal or even stationary - 
lined up along position angles that vary with 
time. These cover a range over more than 90 
(|Tatevama & Kinghaml 12004 lAgudo et allboiol) . 
which is amplified by projection effects owing to 
the small angle of the compact jet to the line 



of sight, 3?2 ± 0?9 IJorstad et all (|2005l ). The 
mean over a 7.5-year period was —107° (mea- 
sured from north, with —90° defined as due west) 
(jTatevama &: Kinghamll2004l ). while the mean po- 
sition angle of the X-ray jet within 8" of the nu- 
cleus is —109°. This implies that, independent of 
the direction of motion on parsec scales, the flow 
bends to collimate along a well-defined direction 
on kilopars ec scales, as has bee n observed directly 
in 3C 279 (Irloman et al.ll2003l) . 



3.2. 



Upper Limit to Unbeamed Radio Lu- 
minosity 



The parsec-scale superluminal m otion, with ap- 



paren t velocities as high as 18c (jJorstad et al 



20051 ). indicates that bulk motions in the jet are 
highly relativistic inside the nucleus. IJorstad et al. 
(|2005l ) combined the apparent speeds with the 
rate of decline of flux of moving knots to derive 
a Doppler factor <5 = 18.9 ±6.4, a bulk Lorentz 
factor r = 16.5 ± 4.0, and an angle between the 
axis and line of sight of 3?2 ± 0°9 for the compact 
jet. The one-sidedness of the extended X-ray and 
radio emission implies that relativistic bulk mo- 
tion continues at least out to an angular distance 
of 20" from the nucleus. 

We cannot determine whether any of the radio 
emission that we detect is unbeamed. However, 
given the limit to the sensitivity of our radio im- 
ages, it is possible that the two outermost features 
seen in Figure [5] are located in a very faint lobe 
or other region where there is only non-relativistic 
bulk flow. We measure a flux density at 1.425 GHz 
from this region to be 4.6 ± 0.9 mJy, which gives 
a 2-cr upper limit to the unbeamed extended radio 
power at 1.4 GHz of < 2xl0 24 WHz" 1 . We there- 
fore confirm that this BL Lac object is within the 
luminosity range expec ted for FR 1 radio sources 
dFanaroff fc Rilevlll974h . 



3.3. Spectral Energy Distribution 

Figure|8]presents the multifrequency continuum 
spectrum of the extended jet. The measurements 
include the entire length of the detected extended 
X-ray jet, i.e., from 3" to 20" from the nucleus 
in order to provide a signal-to-noise ratio at 4.7 
and 14.9 GHz sufficiently high for analysis and to 
allow comparison with the low-resolution IR im- 
ages. The figure reveals that the X-ray emission 
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lies slightly above an extrapolation of the radio 
spectrum, while the 3.6 /im IR and 600 nm opti- 
cal upper limits fall well below the extrapolation. 
In fact, the 600 nm upper limit is comparable to 
the X-ray flux measurement. This is contrary to 
models in which the X-rays arise from synchrotron 
radiation from electrons with energies in the TeV 
range . In such models dDermer k, Atovan _[200S 



20041 : IJester et all 120061 lUchivama et all 
the value of vF v (X-ray) can exceed z/F„(opt), but 
the flux densities cannot unless the electron energy 
spectrum is extraordinarily shallow, Nfa) oc 7~ s 
with s w 1, where 7 is the electron energy in units 
of m c 2 . In this case, essentially all of the electron 
energy would be confined to TeV energies. Fur- 
thermore, the spectral index in the X-ray region 
should then be a(X — ray) ~ 0, while our mea- 
sured value is significantly higher. 

Having eliminated the synchrotron model, we 
construct the SED for the alternative model, 
IC/CMB, by integrating over both the electron 
energy distribution [assumed to be of the form 
AT (7) = 7V 7~ (2q+1) ] and the spectrum of the 
CMB to calculate the X-ray flux density as a 
function of photon energy. In order to avoid vio- 
lating the upper limit at 600 nm, we cut off the 
electron energy distribution below an energy 7 m ; n . 
We present two SEDs for this model in Figured! 
representing the lowest and highest allowed values 
of the product S^ m i n that provide good fits to the 
observed SED. 

4. Discussion 

4.1. Physical Parameters of the Extended 
Jet 

We use the formu l as der ived in the Appendix of 
Jorstad &; Marscherl (|2004r ) to determine the phys- 
ical parameters of the extended jet of O J287 under 
the assumptions that the X-ray emission is from 
the IC/CMB process and that the magnetic and 
particle energy densities are equal, with random 
kinetic energies of the protons and electrons each 
contributing 25% of the total energy density. We 
do this for knot J4, 8" from the nucleus, which is 
well resolved from the other features in both the 
X-ray and radio images. The knot is nearly cir- 
cular (see Fig. [IJ, with a mean extent of a = 1"A. 
We assume that the radio and X-ray spectral in- 
dices are the same, and carry out the calculations 



for the value that is most consistent with the data 
at both wavebands, a = 0.7. 

The above values of a and a, plus the flux densi- 
ties listed in Tabled] yield a Doppler factor (5 = 8.7 
and magnetic field B = 5.4 fiG for the case 7 m i n = 
7 that represents the leftmost curve in Figure [U 
while for the rightmost curve S = 7.9, B — 4.9 /iG, 
and 7 m in = 40. In order to derive the minimum 
total jet power, we assume equipartition between 
the magnetic and electron energy densities. We 
then multiply the magnetic energy density B 2 /8ir 
by 4-7T times the square of the luminosity distance, 
the cross-sectional area 7r(a/2) 2 and speed (« c) of 
the jet, and the square of the bulk Lorentz factor 
(with the approximation that r ~ 5) to transform 
to the host galaxy rest frame. We double the re- 
sult, since the magnetic field supplies only half of 
the total energy under the equipartition assump- 
tion, to arrive at a total power of 1-2 xlO 45 ergs 
s _1 . This is similar to that derived from a de- 
tailed analysis of very long baseline interferometric 
images for the par sec-scale jet of the FR 1 radio 
galaxy 3C 120 by iMarscher etaD (|2007h . Both 
analyses, however, ignore the rest mass of the pro- 
tons, which multiplies the derived kinetic power by 
~ 1800/7~ in , where / is the fraction of the posi- 
tively charged species that are protons rather than 
positrons. The total kinetic power is therefore 
- 1 x 10 45 (l+50/) ergs s" 1 for the value 7 min = 40 
that minimizes the power requirement. This cor- 
responds to ~ 3Mg^ g times the Eddington lumi- 
nosity of a black hole with mass 1 x 10 8 AfBH,8 -Wo- 
lf the mass of the black hole that powers OJ287 is 
as high as 1.8 x 10 10 Mp, , as in the binary model 



of IValtonen et al.l (|2008[ K then the total kinetic 
power is much less than the Eddington value. The 
jet power can be sub-Eddington even if the black- 
hole mass is ~ 10 9 Mq. We therefore find that the 
EC/CMB model for the X-ray emission produces 
reasonable physical parameters, especially if the 
value of 7 m in near the maximum of ~ 40 allowed 
by the SED. 

4.2. Bending, Length, and Deceleration of 
the Jet 

We note that the derived value of the Doppler 
factor in knot it J4 8" from the nucleus is a factor 
of 2 -3 lower than that found in the parsec-scale 
jet (j Jorstad et al.l 120051 ). This implies that the 
jet flow either decelerates on scales of hundreds 
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of kiloparsecs or varies over time-scales of thou- 
sands of years. The main bending of the jet ap- 
pears to occur beyond 8" from the nucleus. If the 
jet remains as straight in three dimensions as it 
appears projected on the sky, as is the most likely 
case, then we can adopt the same value for the an- 
gle of the axis to the line of s ight as that found o n 
parsec scales, 6 = 3?2 ± 0?9 (fjorstad et alJl2005l ). 
The deprojected length of the straight portion of 
the jet is then 6401^40 kpc. In order for a jet at 
an initially small viewing angle to appear to bend 
by more than that angle, it must curve in three 
dimensions by at least the same angle as 9. Such 
bending should require a shock at the boundary 
of the jet opposite the side t oward which it curves 
( Courant fe Friedrichs Il948h . Xnot J3 could rep- 
resent such a shock. If the total curvature of the 
jet in three dimensions is near the minimum re- 
quired to produce the apparent bending by 55°, 
A9 > 3?2, then the final viewing angle increases 
to ^7°. With these values, we estimate the total 
deprojected length of the X-ray jet to be > 1 Mpc, 
with the radio jet extending even farther. 

The Doppler factor that we have obtained from 
our analysis is still quite high, indicating that the 
jet flow remains highly relativistic out to hundreds 
of kiloparsecs from the nucleus. This is contrary 
to the deceleration to sub-relativistic speeds and 
disruption on scales of tens of kiloparsecs that has 
been inferred to be typical of F R 1 radio galaxies 
l|Lainelll99d lLaing et al.lll999h . However, OJ287 
cannot be typical, given the highly relativistic na- 
ture of the jet on parsec scales, with T = 16.5±4.0, 
compared wi th a more common valu e closer to 
r ~ 1 (see iLister fc Marscherl Il997t l. The ex- 
tension of the jet beyond 1 Mpc therefore agrees 
with the theoretical expectation (see §1) that more 
highly relativistic (and presumably more highly 
supersonic) flows can penetrate the interstellar 
and intergalactic medium out to greater distances 
before being slowed down and disrupted. 

4.3. Knotty Structure of the Jet and 
Shock Models 

The bulk of the extended X-ray emission seen in 
Figure [1] emanates from several bright knots. Fig- 
ure [7] shows that the knot to inter-knot flux ratio is 
as high as ~ 10:1 in X-rays and even greater at 4.7 
GHz. A successful emission model must be con- 
sistent with this morphology. In fact, such knotty 



structure has been used as an argument against 
the IC / CMB model on the grounds that the ener- 
gies of the electrons that scatter the CMB photons 
to the X-ray range are low, 7 ~ 40-130 in the case 
of OJ287. Because of this, the electrons have life- 
times to radiative energy losses that are too long 
to confine the emissio n to knots of length 5-10 kpc 
(Sta warz et al. I l2004t) . This is true even for elec- 
trons radiating at centimeter wavelengths, whose 
Lorentz factors are quite high, 7 ~ 10 4 for mag- 
netic fields in the several micro-Gauss range. The 
radiative- loss lifetimes (mainly from IC / CMB ra- 
diation) of these electrons are > 10 6 yr. While gra- 
dients in magnetic field can cause synchrotron in- 
tensity contrasts, the knot to inter-knot flux ratio 
is too high for this to be the sole cause. The knotty 
structure must therefore involve another process. 
We now consider hydrodynamical effects that can 
potentially produce emission features with sizes 
< 10 kpc. 

4-3.1. Standing Shock Model 

Hyd rodynamical effects that oc cur in the flows 
of jets (jCourant fc Friedrichsl[l948T) can cause spa- 
tial variations in both density and bulk veloc- 
ity. Interactions with the external medium form 
shocks oblique to the jet axis t hat decelerate, com - 



press, and bend the jet flow (jGomez et alJ ll997). 



as well as energize relativistic particles. In be- 
tween the shocks are rarefactions, which increase 
the Lorentz factor T at the expense of lower en- 
ergy density u c tx n 7 m i n of the electrons, where 
n c is the number density of electrons. On the 
other hand the energy density of relativistic elec- 
trons is inversely related to T; e.g., in the case 
of adiabatic compression and rarefactions and an 



ultra-relativistic equation of state, V 



-1/4 



(l-2o)-|_i/ 4 



Because the IC/CMB intensity 



-1/4 



[N % 

depends on these factors as Iic/cmb °c N S 2 T 2 
NqT 4 for a small viewing angle, the r a JV ( 
relationship would limit the knot to inter-knot in- 
tensity ratio to a value near unity if compression 
were the only physical effect. We therefore need 
to include non-adiabatic effects of shocks to ex- 
plain the presence of distinct knots of IC/CMB 
emission. 



We use t he formulas given inlLind fc Blandford 
(|l985l ) and ICawthorne fc Cobbl (jl990h to calcu- 
late r and the angle of deflection £ of the flow 
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velocity behind the shock, as well as the ratio of 
post-shock to pre-shock densities. These depend 
on the upstream value of T and the angle of in- 
clination 77 of the shock front relative to the up- 
stream flow direction. We consider the most ex- 
treme case, in which T up = 16.5, as derived on 
parsec scales, is adopted for the upstream flow 
and Tdown ~ 8 downstream of the shock front, 
as in knot J4- This occurs for 77 = 7° and causes 
a deflection £ = 4° that is adequate for a large 
bend to appear in projection, as observed. We 
note that, since it is the velocity component that 
is transverse to the shock front that decreases, 
the shock must be oriented in the direction of 
the bend, which can be seen to be the case for 
knot J3 in Figure [TJ We calculate the upstream 
to downstream density contrast to be a factor of 5 
in this model. If we include compressional heating 

1 /3 

by the shock, so that 7 m i n oc n c , we find that 
A*o oc ^cTm'in 01 2 3 increases by a factor 
of 11 from the compression for our adopted value 
a = 0.7. However, the decrease in T across the 
shock offsets this when calculating the expected 
IC/CMB intensity, since /ic/cmb x NoT 4 . In or- 
der to explain the intensity contrast, the shock 
would need to heat the electrons more strongly 
by another proce ss, probably diffusive particle 
acceleration (e.g., Ellison. Baring, fc Jones 19961 : 



iMeli fc Quenbvll2003l) . In order to bring the knot 
to inter-knot intensity ratio to the observed value 
of ~ 10 by increasing Nq, the minimum (and there- 
fore the mean) relativistic electron energy would 
need to increase by a factor of ~ 16, i.e., 8 times 
more than by compression alone. The required in- 
crease in 7 m i n could be reduced by a factor of 1.6 if 
T up were 10 (roughly the minimum Lorentz factor 
that is compatible with a shock having r<jown = 8) 
instead of the parsec-scale value of 16.5. 

Since bending and deceleration of the jet near 
the X-ray knots can be caused by interactions with 
the external medium that set up standing oblique 
shocks that deflect the flow, a shock model can po- 
tentially explain both features of the jet in OJ287. 
This conclusion is supported by the crude 1.425 
GHz polarization map displayed in Figure |H1 (Un- 
fortunately, the noise level of the polarized inten- 
sity is too high even to estimate the position an- 
gle of the polarization after statistical bias is taken 
into account.) The highest polarization is between 
knot J3 and the southeastern side of J4, where the 



jet first bends noticeably, as well as in J5. The 
polarization is ~ 10% in Knot J2. This conforms 
with the expectation for a shock, which partially 
aligns a chaotic ambient magnetic field along the 
direction of the shock front. 

4.3.2. Moving Shock Model 

Moving shocks, caused by major disturbances 
in the injection at the central engine, are a vi- 
able alternative to stationary ones, since the elec- 
tron energies and densities, as well as the bulk 
Lorentz factor, all increase behind the shock front 
relative to those of the undisturbed flow. Such 
shocks are bounded by rarefactions on the side to- 
ward the nucleus, since the shocked flow is faster 
than the undisturbed flow that was later injected 
into the jet. We note that this scenario is in 
fact similar to the intermittent jet proposal by 
Stawarz et al.l (|2004f ). since it involves a major en- 
hancement in the jet flow velocity and/or injected 
energy density that persists for thousands of years. 
Although the direction of elongation of the knots 
seen in Figure [1] is not always transverse to the 
flow, we will simplify by approximating that the 
shock front is perpendicular to the velocity vec- 
tor. In the rest frame of the shock front, marked 
by a prime, the ratio of density ahead of (sub- 
script "2") to that behind ("1") the front equals 
n 2 /m = (nffiVf r'^), where is the velocity in 
light units (e.g., ISokolov. Marscher. fc McHardv 
2004). We can obtain a consistent solution if, 
in the (unprimed) rest frame of the host galaxy 
of OJ287, Ti = 5 for the pre-shocked ambient 
jet and r s = 11 for the shock front. Relativis- 
tic addition of the velocities gives r 2 = 7.9 for 
the shocked plasma if /3 2 = 1/3, as appr opriate 
for a n ultra-relativistic equation of state (jKonigl 
1980f ). Then n 2 /n 1 = 2.5 and Fi/T 2 = 1.6, which 
is also approximately <5i/c>2. If the electrons are 
heated only by compression, the IC/CMB inten- 
sity, Zic/gmb oc NqT 4 oc n G 7^ n r 4 , behind the 
shock is still ~ 20 times the ambient value. 

Either a standing or moving shock model can 
therefore account for the knotty structure, al- 
though an oblique standing shock is favored, since 
it also explains the deceleration and curvature of 
the jet. Within the model, the slight displacement 
of the radio and X-ray intensity peaks (Fig. [7]) im- 
plies that either the magnetic field or the density 
of electrons with 7 ~ 10 4 is higher near the shock 
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front than 30-50 kpc downstream. In the mov- 
ing shock model, the region of enhancement in B 
and/or No would need to be 30-50 kpc closer to 
the nucleus than the location of the shock front. 
This could be caused by turbulence that develops 
and amplifies the m agnetic field behind the shock 
( Mizuno et al. 2010t ) and possibly also accelerates 
particles. The displacement of the various knots 
from the jet axis implies that the shocks cover only 
a fraction of the jet cross-section. 

4.4. Comment on Spine-Sheath Models 

The strong, one-sided extended X-ray jet in 
OJ287 is best explained if the flow velocity re- 
mains highly relativistic at least out to hundreds 
of kiloparsecs from the nucleus. We detect bright 
knots and ambient X-ray emission across ~ 2" 
(9 kpc) on both sides of the centroid of the jet 
along much of its length. The length-to width ra- 
tio within 8" of the nucleus is only ~ 2, which 
becomes ~ 30 when de-projected, corresponding 
to a full opening angle of ~ 2° when deprojected. 
(We note that this is similar to the open- 
ing ajjgh3_of_yie_£arsec-scale jet derived 
by Ijorstad et al. 20051 . .) Therefore, either the 
spine is not an ultra-narrow pencil beam or its di- 
rection varies, as suggested by the placements of 
the knots. Furthermore, the similar opening angle 
on parsec and kiloparsec scales supports our con- 
clusion that the extended jet remains relativistic, 
since FR 1 jets that decele rate tend to broa den 
considerably on such scales ( Laing et al. 19991 ). 



5. Conclusions 

Our Chandra observations have revealed quite 
strong X-ray emission from the extended jet of the 
BL Lac object OJ287, which has the unbcamcd 
radio luminosity of an FR 1 source. The X-ray 
emission persists out to hundreds of kiloparsecs 
from the nucleus, while weak radio emission can 
be seen beyond 1 Mpc. This demonstrates that 
the main requirement for strong extended X-ray 
emission and an extremely long jet is a highly rel- 
ativistic outflow, independent of the class of radio 
source. 

The observed SED of the extended jet of O J287 
eliminates synchrotron emission as the source of 
the X-rays unless the electron energy distribution 
differs from those in other blazars by concentrating 



essentially all of the energy within the TeV range. 
The highly relativistic flow speed observed in the 
parsec-scale jet and the relativistic beaming im- 
plied by the one-sidedness of the kiloparsec-scale 
jet support an IC/CMB model for the extended 
X-ray emission. Since the relatively low-energy 
electrons involved have extremely long radiative 
lifetimes, the knotty structure of the jet requires 
shocks that compress the plasma and rarefactions 
that lower the density again. The model that ex- 
plains the widest range of observed properties in- 
volves a standing shock at an angle of 7° to the 
upstream flow, which can explain the bending be- 
yond 8" as well as the knot to inter-knot intensity 
ratio and decrease in bulk Lorcntz factor. The 
knotty structure can also arise from variations in 
the jet flow on time-scales of thousands of years. 
This results in moving shock waves followed by 
rarefactions. The knots lie at various positions 
relative to the jet axis, which implies that, rather 
than confinement to an extremely narrow spine, 
the structure subtends a more normal (for rela- 
tivistic jets) full opening angle of ~ 2°, similar 
that that found on parsec scales. 

We selected OJ287 because of its extended, al- 
though very weak, radio jet and its high appar- 
ent superluminal motion on parsec scales. Ex- 
ploration with Chandra of the extended X-ray 
emission around other BL Lac objects with very 
weak kiloparsec-scale radio jets detectable with 
the EVLA could discover more structures with 
high X-ray to radio flux ratios. Observations of 
such sources would further our knowledge of the 
physical development of relativistic jets as they 
propagate through and beyond their host galax- 
ies. 
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Table 1 
Observations 



Telescope 


Instrument 


Band 


UTstart 


UT £nd 


Exposure 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Chandra 


ACIS-S 


0.2-10 keV 


2007-12-27 17:02:36 


2007-12-28 07:15:20 


50.4 ks 


HST 


WFPC2 


600.1 nm 


2007-12-18 21:43:39 


2007-12-19 01:03:39 


2760 s 


Spitzer 


IRAC 


3.6 


2007-11-25 23:30:45 


2007-11-26 00:26:41 


2592 s 


Spitzer 


IRAC 


5.8 fim 


2007-11-25 23:30:45 


2007-11-26 00:26:41 


2592 s 


VLA 


A- array 


1.425 GHz 


1993-01-04 05:30:00 


1993-01-04 14:30:00 


170 m 


VLA 


A- array 


4.710 GHz 


2003-08-31 14:00:00 


2003-08-31 23:00:00 


200 m 


VLA 


B-array 


14.940 GHz 


2007-11-08 12:00:00 


2007-11-08 15:00:00 


150 m 



TABLE 2 

X-Ray and Radio Measurements 



Knot 


$ (X-ray) 1 


F„(1.4 GHz) 2 


F„(4.8 GHz) 2 


F„(15 GHz) 2 


Dimensions 3 


o(X-ray) 4 o(radio) 4 


Jl 




2.0 ±0.2 






< 0.5 




.12 


6.95 ±0.40 


2.2 ±0.3 






1.5 x 1.2 


0.44 ±0.21 


J3 


3.40 ± 0.28 


< 0.15 






1.3 x 1.3 


0.56 ±0.40 


J4 


6.16 ±0.38 


2.9 ±0.2 






1.3 x 1.5 


0.72 ±0.21 


J5 


7.06 ± 0.44 


1.6 ±0.4 








0.68 ±0.22 


J1-J5 5 


30.6 ±0.88 


10.9 ±0.5 


4.13 ±0.43 


1.43 ± 1.03 


1.3 


0.61±0.06 0.8±0.1 


Outer 6 




4.6 ±0.9 











1 Photon counts from 0.2 to 6 keV; units: 10 — 3 ph s — 1 ; no value for Jl because it is confused with the PSF of the nucleus 
2 Flux density in mjy; signal-to-noise ratio is insufficient to determine for individual knots at 4.8 and 15 GHz 
3 Full width across feature at half intensity parallel X perpendicular to axis, in arcseconds 

4 Spectral index, defined such that flux density F v oc v~ Q ; signal-to-noise ratio is insufficient to determine a for some knots, 
especially at the radio frequencies 

5 Includes emission between knots 

6 Region beyond detected X-ray emission out to last radio feature seen in Figure [5] only apparent at 1.4 GHz 
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Fig. 1 .— CXO ACIS-S3 X-ray image of OJ287 
at 0.2-6 keV convolved with a Gaussian kernel of 
FWHM= 0"5. The X-ray knots are labeled. 
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Fig. 2.— HST optical image of OJ287 (gray scale) 
after PSF subtraction, with X-ray (yellow) con- 
tours superposed. Contours are in factors of 2 
starting at 2.075 erg s -1 cm -2 beam -1 , where the 
convolving beam is the same as for Figure [T] The 
image is rotated by ~ 17° clockwise from that of 
the images in the other figures. 




Fig. 3.— SST image at 3.6/im of OJ287 (gray 
scale) after PSF subtraction, with X-ray (yellow) 
contours superposed. Contour levels are as in Fig- 
ure El 
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Fig. 4.— VLA images of OJ287. The global to- 
tal intensity peak is 2.87 Jy/beam at 14.94 GHz. 
The images are convolved with a circular Gaussian 
beam of dimensions lxl arcsec 2 . Contours are in 
factors of 2 starting at 0.01% of the global peak in- 
tensity of 2.87 Jy beam -1 (reached at 14.94 GHz). 
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Fig. 7. — Intensity profiles along the paths marked 
in Figure O Left: X-ray intensity (solid curves) 
and radio intensity at 4.7 GHz (dashed curves) 
along the jet axis. Right: X-ray intensity trans- 
verse to the axis, proceeding from south to north. 



Fig. 5. — Superposed images of OJ287, smoothed 
to enhance the large-scale structure of the jet: 
X-ray (0.2-6 keV, false-color, smoothed with a 
FWHM= 1".5 circular Gaussian function) and 
radio (21 cm, contours, with circular Gaussian 
restoring beam of FWHM= 1'.'5 approximately 
corresponding to natural weighting of the VLA 
fringe- visibility data) . 
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Fig. 6. — Contour map of X-ray intensity, with 
curves and lines indicating the paths correspond- 
ing to the intensity profiles of Figure [7] The X-ray 
image is smoothed with a circular Gaussian restor- 
ing beam with FWHM size 0'.'5, while the contour 
levels are as in Figure [2] 
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Fig. 8. — Continuum spectrum of OJ287. The 
X-ray data are displayed by the flux density and 
uncertainty at 2.4 x 10 17 Hz (1 keV), with solid 
lines corresponding to F v H-a at other frequen- 
cies. The dotted line gives the extrapolation of the 
radio spectrum. The dashed lines correspond to 
the spectrum expected for inverse Compton scat- 
tering of CMB photons for the two extreme values 
of the Doppler factor times the minimum electron 
energy that fit the X-ray and optical data. 
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Fig. 9. — Map of degree of linear polarization of 
OJ287 at 1.425 GHz. The superposed X-ray con- 
tours are the same as in Figure [2j The scale on 
the bottom is in percent. 



